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I.  INTRODUCTION 

In  a  recent  paper,  the  multiple  scale  Fourier  transformation  has  been  applied  to  describe  the 
far-field  asymptotic  behaviour  of  nonlinear  dispersive  waves  (Jeffrey  and  Kawahara  [1]).  It 
was  shown  that  the  use  of  the  multiple  scale  Fourier  transformation  can  systematize  a  nonlinear 
asymptotic  perturbation  analysis,  and  that  the  perturbation  analysis  can  then  be  reduced  to 
simple  manipulations  with  respect  to  Dirac  delta  functions  involving  the  multiple  scale  wave- 
number  and  frequency  space. 

In  the  present  paper,  we  apply  the  multiple  scale  Fourier  transformation  to  oscillation  and 
wave  problems  which  include  damping  or  dissipative  effects.  When  these  effects  are  included  in 
a  system,  the  solution  in  the  physical  space  may  involve  exponential  changes,  but  exponential 
functions  do  not  allow  a  Fourier  transformation  in  the  ordinary  sense.  However,  we  can  intro¬ 
duce  a  Fourier  transformation  for  exponential  functions  if  we  extend  the  notion  of  a  Fourier 
transformation  to  include  generalized  functions.  A  Fourier  transformation  of  an  exponential 
function  can  then  be  represented  by  means  of  a  delta  function  in  a  complex  Fourier  space 
(Bremermann  [2],  Challifour  [3]).  Therefore,  in  terms  of  an  appropriate  definition  of  the  delta 
function  in  a  complex  Fourier  space,  we  can  proceed  with  a  perturbation  analysis  by  means  of 
simple  manipulations  with  respect  to  delta  functions,  even  though  the  system  under  consideration 
may  include  damping  or  dissipative  effects. 

As  a  simple  illustrative  example,  we  consider  the  perturbation  approximation  of  a  linear  oscilla-  ;v 
tion  problem  involving  small  damping.  The  essential  idea  can  be  fully  explained  in  terms  of 
this  example,  and  thereafter  applications  to  general  problems  will  be  a  straightforward  matter. 

In  Section  2,  we  start  by  solving  the  problem  in  terms  of  the  derivative  expansion  method  in 
order  to  make  clear  the  correspondence  between  the  physical  space  representation  and  the 


*  Permanent  address:  Department  of  Engineering  Mathematics,  The  University  of  Newcastle  upon  Tyne,  NE1  7RU, 
England. 
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Fourier  space  representation.  In  Section  3,  we  show  how  the  Fourier  transformation  should  be 
extended  so  as  to  include  exponential  changes.  Once  the  extension  is  introduced,  the  calculations 
become  straightforward,  as  is  shown  in  the  latter  half  of  Section  3.  Some  final  remarks  are  made 
in  Section  4,  and  an  introduction  to  the  notion  and  use  of  a  delta  function  in  a  complex  space  is 
given  in  the  Appendices. 

2.  SOLUTION  BY  MEANS  OF  THE  DERIVATIVE  EXPANSION  METHOD 


We  consider  the  following  linear  oscillation  problem  involving  small  damping 

d2/  ^  ,  d/- 
d 7  +  2edt  +  f  =  °’ 


where  the  initial  conditions  are  given  by 

m  =  a, 


d/(0) 
d  t 


=  0. 


(2.1) 


(2.2) 


This  problem  admits  the  exact  solution 

=  ae_t’[cos(l  —  e2)1/2f  +  e(l  —  £2)_1/2sin(l  —  e2)1/2r].  (2.3) 

First  of  all,  let  us  solve  this  problem  in  terms  of  the  derivative  expansion  method.  Introducing 
the  multiple  time  scales  tm  =  emt  (m  =  0, 1 , . . . ,  M)  and  the  expansions 


/U;e)  =  Z 


H  =  0 


d  _  "  _  d 
*  ~  m?o  K 

into  (2.1)  and  (2.2),  we  find  the  set  of  perturbation  equations 

S% 


K 


+  fo  =  0. 


eJk+f  +2-^_h-2^  =  0 

+  Ji+2dt0dt,ldtn  u’ 


8%  ,  f  ,  ~  ~  d2fi  ,  8% 


K 


+  f2  +2 


_l_  2  , 

8t08tt  8t08t2  '  8tf 


and  the  boundary  conditions 

/o(0 . 0)  =  a,  /,( 0, ..  .,0)  =  f2(  0, ...  ,0)  =  0, 

5/o(0,...,0) 


=  Q  3/1(0.....0)  a/o(0,...,0)  _  Q 

ot0  8tt 

cJ2(  0.....0)  a/,(0,...,0)  cf0( 0.....0) 

a.  +  a,  a,  —  U,  . .  . 


dt2 


(2.4) 

(2.5) 

(2.6a) 

(2.6b) 


+  2fi  +  2ra  =  0’  ••• 

ot0  8tt 


(2.6c) 

(2.7a) 

(2.7b) 
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The  solution  of  (2.6a)  is  given  by 

f0  =  A0  e'*0  +  c.c.,  (2.8) 

where  A0  is  a  complex  function  of  the  slow  variables  f,,  t2, . . . ,  and  c.c.  denotes  the  complex 
conjugate  of  the  preceding  term. 

It  follows  from  (2.6b)  that 


dt20 


izr  +  A  =  ~2i(~^  +  Ao)e‘'°  +  cc- 


The  non-secularity  condition  for  (2.9)  gives 

and  the  non-secular  solution  is  then  given  by 

fi  =  Ai  e"0  -I-  c.c. 

Proceeding  to  (2.6c)  we  obtain  the  non-secularity  condition 


/ £A.  \  f~.<'A0  c2A0  r/40\  n 

2i(  — -1  +  A.  -I-  [2i ^-2  +  - r- j2  +  2-rr  =  0. 

Vat,  v  V  dt\  dtJ 


Solving  (2.10)  we  get 

A0(tvt2,...)  =  B0(t2,...)e-'‘, 
and  substituting  this  into  (2.12),  we  obtain 


+  A  =  Jdlo  +  LBV-,l 
at,  1  \dt2  2  °)e  • 


(2.9) 

(2.10) 

(2.11) 

(2.12) 

(2.13) 

(2.14) 


The  solution  At  for  (2.14)  becomes  bounded  (or  /,//„  -*  0  as  t  -»  oo)  if 


i  ~  0 

(2.15) 

and  is  then  given  by 

/t,(fp  f2,...)  =  B,(t2,...)e‘'‘. 

(2.16) 

From  (2.15)  we  obtain 

(2.17a) 

so  that 

v40(f2,r2 _ )  =  e0(t3 — Je-'*-"-2. 

(2.17b) 

The  next  higher  order  condition  for  the  non-secularity  is  found  to  be 

JdA2  A  \  /.CA.  d2A.  ?A.\  J.cA0 

2fl  +  A2  )  +  (2i  — ^ 1  +  — r-j1  +  2-x-1 )  +  2  i-,-2  + 

Vat,  V  V  dtJ  V 

£*,♦£)-*•  ai8) 
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Similar  calculations  are  also  possible  for  higher  order  terms,  and  up  to  the  0(e4)-solution  we 
arrive  at 

/  =  E0c“°~,'~i'll2'i'418  +  c.c. 

+  e(Dl  e‘'°~'1_i'2/2  +  c.c.) 

+  e2(C2  e"0-'1  _"2/2  +  c.c.) 

+  e2(B3  e1'0-'1  +  c.c.)  +  e"°  +  c.c.).  (2.19) 

The  boundary  conditions  (2.7)  determine  the  coefficients  of  (2.19)  successively,  giving  rise  to 

E0  =  a/2,  />,  =  -  iaj 2,  C2  -  0, 

B3  =  -  ia/4,  Aa  =  0,  (2.20) 

where  £0,  Dv  C2,  Bv  and  Aa  are  assumed  to  be  constants  with  respect  to  the  slow  variables  equal 
to  or  slower  than  f5,  t4,  r3,  f2,  and  tr  respectively.  Consequently,  up  to  the  O(e4)-approximation 
calculated  here,  we  obtain 


/ 


'^2/2  1/4/8  __)  i  jg  gi/o  ^2/2  g^o  |  g 

2 


.c. 


+  0(£5) 


(2.21) 


=  ?e-»|e"°-„; 

=  a  e_££cos^l  —  4-  e  sin^l  +  ^  sin  t  +  0(£5). 

It  should  be  noticed  here  that  the  solution  (2.21)  is  not  a  simple  expansion  in  e  of  the  exact 
solution  (2.3).  The  expansion  of  the  factor  (1  —  £2)1/2  in  the  exact  solution  reproduces  exactly 
the  terms  in  (2.21)  up  to  the  order  of  approximation  calculated  here.  The  derivative  expansion 
that  avoids  the  secularly  thus  incorporates  partial  sums,  in  the  sense  that  the  perturbation  solu¬ 
tion  so  obtained  is  not  a  simple  power  series  solution  in  e. 

3.  MULTIPLE  SCALE  FOURIER  TRANSFORMATION 

We  shall  now  solve  the  problem  (2.1)  subject  to  (2.2)  in  the  Fourier  space.  The  Fourier  trans¬ 
formation  for  the  function  /  with  multiple  scales  t0,  t tM  is  first  introduced  as  follows: 


(3.1) 


Hoo  C  M  1  M 

</( W0,(0 a>„) exp <|i  £  ayj  da>m, 

-  ao  l  m  =  0  )  m  =  0 

o> wi’  •  •  •  ’  =  (2n)M  + 1  /  J”  *i» •  •  •  > *»#)exP  ~ *  X  0  (3-2) 

At  the  moment  we  assume  that  to0,  to,, . . . ,  are  real. 

In  agreement  with  the  expansions  (2.4)  and  (2.5),  we  expand  the  Fourier  amplitude  as  follows: 


g(u> 


g((o0,  to,, . . . ,  o)J  =  £  engll(to0,(Dl,...,toJ. 


(3.3) 


n  =  0 


It  then  follows  from  (3.1)  that  the  time  derivative  should  be  replaced  in  the  multiple  scale  Fourier 
space  by  the  operation 


i(<W0  +  ECO,  +  e2(o2  +  ...). 


(3.4) 
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The  introduction  of  (3.3)  and  (3.4)  into  (2.1)  and  (2.2)  now  yields 


(“«S  +  1]»0  -  0.  (3-5a) 

(-<ol  +  l)gl  -  2co0(to,  -  i)g0  =  0, 


(~w20  +  l)g2-  2ai0(wl  ~  i)gt  -  (<o\  +  2 to0u)2  -  2 i(ot)g0  =  0,  (3.5c) 


while  the  boundary  conditions  reduce  to 


K 0.  +  W10o)  FI  <K 


m  =  0 

At 

(w0g2  +  i olg1  +  (o2g0)  f]  d«m  =  . . .  =  0. 

m  -0 


(3.6a) 


(3.6b) 


It  will  be  shown  in  the  course  of  the  calculations  that  follow  how  the  necessity  arises  for  the 
extension  of  c o0,  a>,, . . . ,  to  the  complex  plane. 

The  general  solution  to  (3.5a)  may  be  expressed  in  the  form 


0oK)>  °V  •  •  •)  =  _  l)a0(ci>,,a)2, ...)  +  <H<o0  +  l)a*(-«V  -to2,...),  (3.7) 

where  we  have  made  use  of  the  reality  of  f0,  and  *  denotes  a  complex  conjugate.  It  is  readily 
seen  that  the  inverse  Fourier  transformation  of  (3.7)  reproduces  (2.8),  if  we  put 

M  ^  At 

«  I  "m'J  n  d0)m-  (3-8) 

m  =  1  )  m  -  1 

The  substitution  of  (3.7)  into  (3.5b)  leads  directly  to  the  result 

2o)n 

g ,-••)=  2  _  .  t&o  -  OK  -  i)a0{(0v<o2,...) 

cu0  1 

- T3L7<*H)  +  1)K  -  OaS(-<Up  -a>2,...).  (3.9) 

Now  it  may  be  seen  from  the  expression  (3.9)  that  gx  becomes  singular  at  <u0  =  ±1,  because  the 
denominator  to*  —  1  becomes  zero  at  the  origin  of  the  respective  delta  functions  <5(co0  -  1) 
and  <5(co0  +  1). 

These  singularities  can  be  eliminated  if  we  can  set 


A>^1*  ^2’  ■  ■  •) 


-  j; 


a0(w,,o>2,...)cxp 


S{o)0  -  1) (to,  -  i)a0(o)t,o)2, . . .)  =  0, 

S(oi0  +  l)(co,  -  Oa*(-to,,  -to2,...)  =  0, 


at  to*  =  ±  1,  or  if  we  can  set 

a0(<ov  to2, . . .)  -  <5(to,  -  i)  b0(a>2, . . .),  3  j 

a*(-to,,  -to2,...)  =  S(wi  -  i) b*(— to2, . . .). 

So  for  to0,  a>,, . . . ,  have  been  assumed  to  be  real,  so  that  to,  is  real,  and  consequently  the  ex¬ 
pression  <5(co,  -  i)  does  not  possess  a  well-defined  meaning.  As  is  easily  seen  from  (2.10)  and 
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(2.13),  the  expression  d(c o,  —  i)  corresponds  to  the  exponential  change  e-'1  in  the  physical  space. 
The  exponential  function  is  not  absolutely  integrable,  and  so  it  has  no  Fourier  transformation 
in  the  ordinary  sense.  However,  we  can  apply  the  Fourier  transformation  for  generalized  func¬ 
tions  to  the  present  example.  It  is  then  necessary  to  extend  the  Fourier  space  w0,  a>,, . . . ,  into  a 
complex  space. 

It  is  possible  to  introduce  the  Fourier  transformation  of  an  exponential  function  by  appeal  to 
the  theory  of  Fourier  transformations  for  generalized  functions  (Challifour  [3]).  The  following 
formulae  that  result  are  useful  for  the  present  analysis: 


^[ew]  =  2nd(a>  +  bi). 

(3.12a) 

^[eial]  =  2nS(o)  +  a). 

(3.12b) 

e*']  =  2nS(o)  +  a  +  biX 

(3.12c) 

where  co  has  been  extended  to  a  complex  space,  and  the  details  of  the  extension  are  shown  in 
Appendix  A.  Formula  (3.12a)  may  be  used  to  define  <5(a>,  —  i)  in  the  present  problem.  The  re¬ 
lation  (3.12b)  indicates  that  a  sinusoidal  change  in  the  physical  space  corresponds  to  a  delta 
function  in  the  Fourier  space,  shifted  by  an  amount  a,  as  already  noted  previously  (Jeffrey  and 
Kawahara  [  1  ]).  Formula  (3.1 2c)  will  be  useful  for  other  general  problems,  as  is  shown  in  Appendix 
B. 

For  practical  calculations,  it  is  convenient  to  assume  merely  that  a>0,  a>,, . . . ,  have  all  been 
extended  to  complex  quantities,  and  thereafter  to  carry  out  calculations  in  terms  of  delta  func¬ 
tions  in  a  complex  space.  For  the  present  problem,  it  suffices  to  introduce  the  complex  quantity 
co,  for  the  delta  function  <5(co,  - i),  but  we  extend  (3.1 )  and  (3.2)  to  make  them  valid  for  all  complex 

"»• 

The  lowest  order  solution  is  then  given  by 

0o(coo,co,, . . .)  =  S(co0  -  1  fflco,  -  i)b0(oj2, . . .)  +  S(oj0  +  l>5(a>,  -  j>*(-co2, . . .).  (3.13) 

The  non-secular  solution  for  g2  that  satisfies 

(«o  “  «,.-••)  =  0,  (3.14) 

is  provided  by 

g ,(0)0,0),,...)  =  6(a>0  -  lja^to,,...)  +  d(o>0  +  . .),  (3.15) 

which  is  equivalent  to  the  solution  (2.1 1)  in  the  physical  space  representation. 

Introducing  (3.13)  and  (3. IS)  into  (3.5c),  we  obtain 

(o>o  -  l)g2(a>0,  Op . . .)  =  2co0(co,  -  i)  [<J(a>0  -  l)o,(a),, . . .)  +  6(cu0  +  l)a*( -<av . . .)] 

+  (to?  +  2o)0o>2  -  2ict>,)[<S(a>0—  lj^to,  -  i)  b0(w2, . . .)  +  6(io0  +  l^fco,  -  0 

x  (3.16) 

Now,  singularities  arise  at  <w0  =  ±  1  when  we  solve  (3.16)  for  g2,  but  they  can  be  eliminated  if 
we  can  set 

2(a),  -  Oa,(a>,, . . .)  +  (1  +  2 o)2)S(a>t  -  i)b(a)2, . . .)  =  0, 

-2(co,  -  i)a*(-cD,,...)  +  (1  -  2 o)2)S((ot  -  i)b*(-(o2,...)  «  0, 


(3.17) 
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This  relation  corresponds  to  (2.12)  in  the  physical  space  representation. 

When  we  solve  (3.17)  for  a,,  a  singularity  arises  at  to,  =  i,  because  of  the  appearance  of  the 
factor  5(a>,  —  /)/( ai,  —  i ).  To  remove  such  singularities  we  can  choose 


(1  +  2a)2)60(a>2, . . .)  =  0,  (1  -  2ai2)fc*(  — oi2, . . .)  =  0, 

(3.18) 

that  is, 

^0^2’ 

oi3, . . .)  =  5(oi2  +  \)c0(ojv  . . .),  b*(  -  (o2,  -  a)3, . . .)  =  5(oi2  -  ^)c*(  -  ai3, . 

..). 

(3.19) 

Then  al(a>l 

, . . .)  is  given  by 

a,(a),, . ..)  =  5(a),  -  06,(a)2, . . .), 
a*(- ai, ,...)  =  5(a),  -  i)b*(  -  o)2, . . .). 

(3.20) 

The  inverse  Fourier  transformations  of  (3.19)  and  (3.20)  reproduce  the  physical  space  representa¬ 
tions  (2.16)  and  (2.17). 

Using  (3.19)  and  (3.1 1)  in  (3.7),  we  find 

y0(a)0,  tov  ...)  =  5(oi0  -  l)5(o),  -  i)S(w2  +  ^c0(co3, . . .)  +  6(w0  +  1)5(01,  -  i)5(a)2  -  ^) 

x  c*(-ai3,...).  (3.21) 

The  inverse  Fourier  transformation  of  (3.21)  gives 

C0(t3, ...)  e"0"“j/2e'*'  +  CJ(t3, . .  .)e-fto+ftl/2  e_,‘,  (3.22) 

where  C0(t3, . . .)  represents  the  inverse  Fourier  transformation  of  c0(a>3, . . .). 

The  next  order  problem  is  given  by 

(  oIq  +  1  )g3  =  2 d(o)0  -  l)[(oi,  -  i)a2  +  5(a),  -  i)  (oi2  +  $>,  +  5(ai,  -  i)  S(co2  +  |)oi3c0] 
— 25(ai0  +  1 )  [(ai,  -  i)a*  +  5(a),  -  i)(a>2  -  ^)6*  +  5(ai,  -  i )5(ai2  +  ^)a)3c*].  (3.23) 
The  non-secular  condition  for  is 

(O),  -  i]a2(o)v . . .)  +  5(0),  -  i)  (ai2  +  ^)&,(a)2, . . .)  +  5(a),  -  i)5(a)2  4-  £)a)3c0(a>3, . . .)  =  0, 

(3.24) 

together  with  a  similar  relation  involving  asterisks.  The  relation  (3.24)  corresponds  to  (2.18) 
in  the  physical  space  representation.  When  we  solve  (3.24)  for  a2,  there  arises  the  singularity 
5(a),  -  i)/(a),  -  i).  It  can  be  removed  if  we  set 

(a)2  +  2)6,(a)2, . . .)  +  5(a>2  +  ^)a)3c0(a)3, . . .)  =  0.  (3.25) 

Similarly,  a  singularity  related  to  5(a>2  +  2)/(ai2  +  ^arises  when  seeking  ft,.  It  also  can  be  avoided 
if  we  set 

a)3c0(a>3,  ••.)  =  (),  (3.26a) 

or.  if  we  set 

c0(ai3,  a>4, . . .)  =  5(ai3)d0(a)4, . . .).  (3.26b) 

We  can  proceed  with  similar  calculations  up  to  any  order  of  approximation  in  a  straightfor¬ 
ward  manner.  For  example,  up  to  the  O(e*)-approximation,  we  obtain 
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y  =  y0  +  £y  i  +  e^2  +  £303  +  «4«4 

=  S((00  -  1)[S((01  -  i)S(co2  +  j)S(a)3)S((04  +  £)<?„  +  &5(cu,  -  i)S(io2  +  $)5((o3)dl 
+  e25((oi  -  i)S((o2  +  j)c2  +  e3d((o1  -  i)b3  +  e4a4]  +  <5(a>0  4-  1) 
x  [S{0)1  -  i)S(w2  -  j]d{a>3}d(a)4  -  ±)e*  +  eS(ws  -  i)S(a)2  -  ^(c o3)d*  +  e2d((o1  -  i ) 

x  d(o)2  —  j)c*  +  e3d (ajj  -  i)b*  +  e*a*].  (3.27) 

The  use  of  the  boundary  conditions  for  this  solution,  together  with  the  assumption  that  e0, 
d„  c2,  b3,  and  a4  are  constants  with  respect  to  the  slow  variables  equal  to  or  slower  than  cu5, 
co4,  co3,  co2,  and  a>r  respectively,  leads  to  a  solution  equivalent  to  (2.21)  in  the  physical  space 
representation.  Thus,  once  the  delta  function  in  a  complex  Fourier  space  has  been  introduced, 
we  can  carry  out  the  perturbation  analysis  systematically  in  the  Fourier  space. 

4.  CONCLUDING  REMARKS 

Sandri  [4]  applied  the  multiple  scale  Fourier  transformation  to  a  problem  involving  expon¬ 
ential  changes.  However,  there  was  an  error  in  his  paper,  since  he  omitted  the  imaginary  unit  in 
his  calculations  on  page  91.  Because  of  this  error,  he  was  unaware  of  the  point  discussed  in  this 
paper  that  the  delta  function  should  be  extended  to  a  complex  space,  corresponding  to  the  Fourier 
transformation  of  an  exponential  function. 

If  the  Fourier  transformation  is  extended  to  a  complex  space  on  a  mathematically  sound  basis, 
we  can  apply  the  multiple  scale  Fourier  transformation  to  perturbation  problems  in  general 
physical  systems,  merely  keeping  in  mind  that  the  multiple  scales  in  the  Fourier  space  are  complex 
quantities.  In  practical  calculations,  the  use  of  the  delta  functions  in  the  complex  space  may 
simplify  the  manipulations  associated  with  the  perturbation  analysis. 

Although  the  example  treated  in  this  paper  is  simple,  the  essential  point  involved  is  well  illus¬ 
trated.  The  extension  presented  in  this  paper  will  benefit  the  Fourier  transformation  method 
when  it  is  applied  to  solve  general  wave  problems  which  involve  dissipative  effects. 
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APPENDIX  A.  FOURIER  TRANSFORMATION  OF  EXPONENTIAL  FUNCTIONS 
Let  us  denote  the  Fourier  transformation  of  /  by 

r  t  f,m)  =  8(0.)  =  J  /(f)  e1"  dt,  (A.1) 
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and  define  its  inverse  by 


g(a))ei0i'da>. 


When  /  and  g  belong  to  the  class  of  L,  or  L2 -functions,  their  Fourier  transforms  exist  in  the  ordinary  sense.  If  this  is 
not  the  case,  we  must  then  introduce  the  notion  of  generalized  functions  together  with  an  extended  definition  for  the 
Fourier  transformation.  Such  generalized  Fourier  transformations  can  be  introduced  in  terms  of  those  for  distributions 
(Bremermann  [2],  Chailifour  [3])  or  those  for  hyper-functions  (Imai  [5]). 

The  distribution  generated  by  /(f)  is  defined  by 

</,4>>  =  J  M<H t)dt,  (A.3) 

where  <f>  is  a  test  function.  Let  T  be  a  distribution,  then  its  Fourier  transformation  is  introduced  as  follows 


when  the  inversion  theorem  becomes 


<^-J(r)),0>  =  -  <T4>>-  (A.5) 

As  an  example,  let  us  consider  that  1  denotes  the  distribution  generated  by  the  unit  function  1.  Then  it  follows  that 

<t>>  =  <(ku,  ei0*>,  -  L  (A.6) 

Thus  we  have  &(6)  =  I.  By  means  of  (A.5)  we  have  *(1)  =  S.  Furthermore,  it  is  easily  shown  that  &  ~l{6)  =  and 
^(1)  =  2nS. 

Let  us  now  consider  the  Fourier  transformation  of  e*  Using  the  formulae  (A.3)  to  (A.5),  we  have 


<^(e*a>),0(a>)>  =  <eia',J^(0,f)> 


Hence  we  have 


=  ^1,  j*  <t>(oj  —  a)e'""  dco^ 

—  <2 7tS,  <p((o  —  a))  =  2n<p{—a) 

—  <2?r^(cu  4-  a),  <£(co)>. 

iHe*-1,  co)  =  2n S{(o  4  a). 


The  formula  (A.8)  has  been  obtained  for  real  a  in  the  above  proof.  We  may  also  give  a  similar  proof  for  a  complex  cu, 
in  order  to  include  the  Fourier  transformation  for  an  exponential  function,  say  eH.  It  is  possible  to  choose  the  test  function 
<t>  such  that  it  vanishes  outside  a  compact  set,  and  thus  the  integral  of  the  product  <f>  titat  with  e*  becomes  bounded  along 
the  appropriate  integration  path  in  the  complex  space.  In  such  a  way,  the  Fourier  transformation  of  e*  may  be  given  as 

•  ,^(e*,a>)  *  2 nd{(o  4  bi ).  (A.9) 

Consequently,  we  arrive  at  the  formulae  (3.12)  given  in  the  text. 

The  Fourier  transformation  of  generalized  functions  can  also  be  introduced  in  terms  of  hyper-function  theory  (Imai 
[5]).  We  do  not  go  into  details  here,  but  notice  only  that  the  formulae  (3.12)  are  also  obtainable  in  terms  of  the  Fourier 
transformation  of  hyper-functions. 


APPENDIX  B.  APPLICATION  OF  THE  DELTA  FUNCTION  IN  A  COMPLEX  SPACE 

Let  us  solve  the  problem  (2.1)  subject  to  (2.2)  in  a  Fourier  space,  assuming  that  e  is  of  order  unity.  We  consider  the 
Fourier  transformation 

/(()=  f  t/Me'-dco,  (B.1) 


f 


V%  ■ 
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where  to  is  extended  to  a  complex  space,  and  y  is  a  suitable  contour. 

Introducing  (B.l)  into  (2.1)  and  (2.2),  we  obtain 

(  —  to2  +  2ie(o  +  l)g{to)  =  0,  (B.2) 

J  g(vj)  dc o  -  a,  J  icogio)  da>  —  0.  (B.3) 

When  we  admit  delta  functions  in  the  complex  cn-space,  the  general  solution  of  (B.2)  can  be  expressed  as  follows: 

g{to)  =  <5(co  -  -  ie)C  4-  <5(o>  +  yT  -  r  -  it)D. 

The  inverse  Fourier  transformation  then  gives 


(B.4) 


At) 


~  j  <5(«j  —  N/1  —  f.z  —  ie)C  clltit  dw 

+  I  <5(cu  +  yj\  -  P  —  ie)D  Qi<ot  dco 

=  e  ~et[C  exp  iyf  1  ~ e 2  t  +  D  exp  -  iyl  -  e2  r]. 

The  boundary  conditions  give  the  relations 

C  +  D  =  n,  iy'T"^"F(C  -  D)  -  e(C  +  D)  =  0, 

which,  in  turn,  give  rise  to 


(B.5) 

(B.6) 


b-K,+JT^7V 


Thus  we  obtain  the  exact  solution  (2.3)  quite  easily  in  terms  of  a  delta  function  in  a  complex  Fourier  space. 
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